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Abstract. Pulsars are amongst the most stable rotators known in the Universe. Over many 
years some millisecond pulsars rival the stability of atomic clocks. Comparing observations of 
many such stable pulsars may allow the first direct detection of gravitational waves, improve the 
Solar System planetary ephemeris and provide a means to study irregularities in terrestrial time 
scales. Here we review the goals and status of current and future pulsar timing array projects. 
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1. Introduction 



The pulsar timing technique (e.g. Lorimer & Kramer 2005 for an overview and Ed- 
wards, Hobbs & Manchester 2006 for a detailed description) has enabled some of the most 
exciting recent results in astrophysics. For instance, the first extra-Solar planets were dis- 
covered orbiting the pulsar PSR B1257+12 (Wolszczan & Frail 1992), the first evidence 
Qh' for gravitational waves was provided by the binary system PSR B1913+16 (Hulse & Tay- 
lor 1974) and the most stringent tests of general relativity in the strong-field regime have 
been carried out using the double pulsar system (Kramer et al. 2006). Even though the 
PSR 131913+16 system has provided evidence for the existence of gravitational waves, 
most researchers regard this as indirect evidence and not a direct detection of such waves. 

Sazhin (1978) and Detweiler (1979) showed that low-frequency gravitational waves will 
induce pulsar timing residual^ that could be detectable with high-precision observations 
of the pulse arrival times. However, with a single pulsar it is difficult, and perhaps im- 
possible, to distinguish between the effects of gravitational waves and many other effects 
such as the irregular rotation of the neutron star, interstellar medium propagation effects 
or inaccuracies in the planetary ephemeris. Hellings & Downs (1983) showed that such 
effects could be distinguished by looking for correlated behaviour in the timing residuals 
y—i ' of many pulsars. In brief, the timing residuals caused by irregular rotation or propagation 
Q\ . effects should be uncorrelated between different pulsars. Irregularities in the terrestrial 
atomic time scale would produce completely correlated timing residuals. The angular 
, correlation function expected for an isotropic, stochastic gravitational wave background 
(Hellings & Down 1983) is shown in Figure [1] (using simulated data described in Hobbs 
et al. 2009). 

The first attempt to create a "pulsar timing array" in which enough pulsars are ob- 
served with sufficient timing precision to search for correlated signals was reported by 
Foster & Backer (1990). More recently Jenet et al. (2005) calculated the number of pul- 
sars, the timing precision and data span required in order to make a significant detection 

f Pulse times of arrival for a given pulsar are compared with a physical model of the pulsar. 
The differences between the measured arrival times and predicted arrival times using the model 
are known as the pulsar timing residuals and denote physical effects that are not included in 
the model. 
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Figure 1. The expected correlation in the timing residuals of pairs of pulsars as a function 
of angular separation for an isotropic gravitational wave background. The solid line gives the 
theoretical curve. The dots are calculated by correlating simulated data sets for 20 pulsars in 
the presence of a gravitational wave background that produces white timing residuals. 

of the expected stochastic background of gravitational waves. In summary, at least 20 mil- 
lisecond pulsars are required, observed for five years with timing precision around 100 ns. 
In this review paper we discuss the current status of various existing pulsar timing array 
projects that aim to achieve this goal. In the future, telescopes such as the Square Kilo- 
metre Array may make gravitational wave detection using pulsars commonplace. The 
current status of this project is described in iJ3] 



2. The pulsar timing array projects 

The main aims of pulsar timing arrays are to search for correlated signals in the pulsar 
timing residuals in order to 1) improve the Solar System ephemeris, 2) look for irregu- 
larities in the terrestrial time standards and 3) make a direct detection of gravitational 
waves. Three projects have been started: 

• The Partes Pulsar Timing Array is the only timing array project (to date) in the 
Southern Hemisphere. This project makes use of the 64-m Parkes radio telescope. Every 
two-three weeks, each of 20 pulsars is observed for approximately one hour using a 
dual-band 10cm-50cm receiver and for another hour using a 20cm receiver. The data 
are processed using the PSRCHIVE suite of software (Hotan, van Straten & Manchester 
2004) that allows polarimetric and flux density calibration. Pulsar timing models and the 
resulting timing residuals are obtained using tempo2 (Hobbs, Edwards & Manchester 
2006). 

• The European Pulsar Timing Array is making precise timing observations of approx- 
imately 15 pulsars. This project makes use of the existing four 100-m class telescopes in 
Europe (Jodrell Bank, Effelsberg, Westerbork and Nancay) and will also include obser- 
vations from the 64-m Sardinian radio telescope after its completion. As well as using 
each of these telescopes as individual instruments, the Large European Array for Pulsars 
(LEAP) project will allow these telescopes to act as a phased array giving the equivalent 
sensitivity of a 200 m telescope. 

• The North American Pulsar Timing Array, known as Nanograv, combines data col- 
lected with the 300-m Arecibo radio telescope and the 100 m GreenBank telescope. These 
telescopes are currently obtaining about monthly observations on 24 pulsars. 
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Figure 2. Simulated timing residuals for PSR J1857+0943 (left panel) and J0437-4715 
induced by an error in the mass of Jupiter of 10 -8 Mq . 

Members of these three projects have agreed to share data to form the International 
Pulsar Timing Array (IPTA) project which hopefully will result in the detection of grav- 
itational waves within 5-10 years. 

2.1. The planetary ephemeris 

The pulsar timing procedures requires that the observed pulse times-of-arrival are cor- 
rected to the Solar System barycentre. This correction is carried out using a planetary 
ephemeris that provides the position of the Earth with respect to the barycentre. Such 
ephemerides are complex and take into account the masses and orbital parameters of the 
planets and over 300 asteroids. If the mass of one planet was slightly incorrect, then the 
position of the barycentre will be incorrectly determined and sinusoidal timing residuals 
will be induced with the period of the planet's orbit (note that a change in the mass 
of Jupiter will affect the entire fit that is carried out to create the planetary ephemeris. 
The change in position of the barycentre is a first order effect and, as the change in 
Jupiter's mass considered here is small, is the only change in the ephemeris that is con- 
sidered. A more detailed study will be published elsewhere.) In Figure [5] we simulate the 
timing residuals of PSR J1857+0943 (left panel) and PSR J0437-4715 (right panel) for 
a change in the mass of Jupiter of 1O _8 M0. The IPTA project is currently combining 
data from a selection of pulsars to place a limit on the error in Jupiter's mass which is 
expected to be more precise than measurements obtained using the Voyager space-craft, 
but not as precise as unpublished Galileo space-craft measurements. In the near future 
pulsar timing arrays should be able to rule out, or detect, the existence of an Earth-mass 
"planct-X" closer than 60 A. U., or a Jupiter-mass object closer than 300 A. U. 

2.2. Terrestrial time standards 

Millisecond pulsars are extremely stable. A measure of stability, a z (Matsakis et al. 1997), 
allows pulsars to be compared with terrestrial time scales. In Figurc[3]we plot the stability 
of PSRs J0437-4715 (crosses) and J1909-3744 (squares). The arrival times for both these 
pulsars were referenced to the Bureau International des Poids et Mesures (BIPM) 2008 
terrestrial time standard. The solid line indicates the stability of the difference between 
the Physikalisch-Technische Bundesanstalt (PTB) and National Institute of Standards 
and Technology (NIST) time standards. The dotted line corresponds to timing residuals 
with a flat spectrum. It is clear that on short time scales the atomic time scales are 
significantly more stable than pulsars, but on time scales of ~ 10 yr and longer the 
pulsars become as stable as the terrestrial clocks. In the near future it is expected that 
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Figure 3. a z (r) for PSR J0437-4715 and PSR J1909-3744 and for a data sequence of atomic 

clock differences: PTB-N1ST. 
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Figure 4. Characteristic strain sensitivity for existing and proposed GW detectors as a function 
of GW frequency. Predicted signal levels from various astrophysical sources are shown (Hobbs 
2008). 



correlated signals in the timing residuals will indicate the presence of irregularities in the 
terrestrial time scales and will allow the formation of a pulsar-based time standard (e.g. 
Rodin 2008). 

2.3. Gravitational wave detection 
The most likely detectable gravitational wave signal is from an isotropic, unpolarised, 
gravitational-wave background caused by binary supermassive black holes at the cen- 
tre of merged galaxies (Sesana et al. 2008). The induced timing residuals induced by 
such a background can be simulated using the TEMP02 software package (Hobbs et 
al. 2009) and lead to upper bounds on the amplitude of any such signal. Jenet et al. 
(2006) gave an upper limit on the energy density per unit logarithmic frequency interval 
of S! 9 [l/(8yr)]ft, 2 ^ 1.9 x 10~ 8 which can be used to constrain the merger rate of su- 
permassive binary black hole systems at high redshift (Wen et al. 2009). A comparison 
between the sensitivity to gravitational waves for different detectors is shown in Figured] 
along with the predicted signal levels from various astrophysical sources. 



3. Future timing array projects 

With current telescopes it is a challenge to achieve a timing precision of ~100ns 
for more than a few pulsars (the best published data set is for PSR J0437— 4715 ob- 
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served for more than 10 yr with rms timing residuals ~ 200 ns; Verbiest et al. 2008). 
Within ~10 yr it is expected that various new telescopes such as the Australian Square 
Kilometre Array Pathfinder (ASKAP) and the South African pathfinder (MEERKAT) 
telescopes will increase the number of known pulsars and provide a larger number of 
telescopes in the Southern Hemisphere for high-precision pulsar timing (Johnston et al. 
2008, http://www.ska.ac.za). Also a new 500m-diameter telescope (FAST) should be 
completed in China (Nan et al. 2006). This telescope will be a highly efficient pulsar 
search and timing instrument and should significantly improve our timing precision for 
a large number of pulsars. 

On a longer timescale, the Square Kilometre Array (SKA) telescope is planned. This 
telescope should be able to observe many hundreds of pulsars with the precision currently 
achieved for only a few pulsars. Such data sets should allow the gravitational-wave back- 
ground and individual sources of such waves to be studied in detail (Cordes et al. 2004, 
Kramer et al. 2004). For instance a high S/N detection of a stochastic background would 
provide a test of general relativity (Lee et al. 2008) and allow a detailed understanding 
of the properties of the sources that form the background. 

4. Conclusion 

Many pulsars are currently being observed as part of pulsar timing array projects. 
Correlated timing residuals allow 1) irregularities in terrestrial time standards to be 
identified, 2) inaccurate planetary masses in the Solar System ephemeris to be corrected 
and 3) the detection of gravitational wave signals. Combining observations from Northern 
and Southern hemisphere telescopes should allow the detection of gravitational waves on 
a time scale of 5-10 years. 
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